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ELECTRONIC  POLARIZABILITIES  AND  STERNHEIMER  SHIELDING  FACTORS 

ABSTRACT 

A  new  method  is  developed  for  determining  the  distortions  (polariza- 
bilities)  induced  in  electronic  distributions  by  valence  electrons  and/or 
crystalline  fields  and  their  effect  (expressed  as  Sternheimer  shielding 
factors)  on  magnetic  and  electric  hyperfine  interactions.  For  illustrative 
purposes  emphasis  is  placed  in  this  paper  on  the  calculation  of  Sternheimer 
anti -shielding  factors  (y^).  Working  within  the  framework  of  the  Hartree- 
Fock  sel f -consistent -  field  formalism,  it  is  shown  that  the  'angular  *  ex¬ 
citations  are  gotten  by  relaxing  the  usual  restriction  that  the  spatial 
part  of  the  one-clectron  functions  be  separable  into  a  radial  function 
times  an  angular  function;  relaxing  the  restriction  that  electrons  of  the 
same  shell  but  differing  in  magnetic  quantum  number  (m*)  have  the  same 
radial  function  yields  the  'radial'  excitations.  To  illustrate  the  method, 
calculations  are  reported  for  several  spherical  ions  (Cl*  and  Cu+)  in  an 
external  field,  but  the  scheme  is  also  applicable  to  the  problem  of  induced 
electric  quadrupole  (and  magnetic  dipole  and  higher  multipole)  distortions 
of  an  ion  by  its  own  aspherical  charge  distribution.  The  problems  of 
orthogonality,  exchange,  and  self-consistency,  which  have  complicated 
applications  of  the  perturbation  method  are  easily  resolved  by  this  approach. 
Further,  since  a  self-consistent  field  procedure  is  followed,  the  distortions 
induced  in  the  inner  closed  shells  by  the  distorted  outer  shells  are  included 
in  a  natural  way  and  by  comparison  with  the  results  of  the  perturbation- 
variation  method  (which  does  not  take  these  into  consideration),  these 
additional  effects  are  shown  to  be  significant. 
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I.  INTRODUCTION 


Electric  quadrupole  and  magnetic  dipole  interactiona  between  atomic 
nuclei  and  outer  electron  distributions  have  been  measured  in  atoms  and 
molecules,  and  in  metals  and  salts  by  a  variety  of  methods  and  most  recently 
by  recoil-less  emission  and  absorption  of  y-rays  (the  MOsabauer  effect). 

As  emphasized  by' Sternheimer 1 ,  the  interpretation  of  these  experiments  (such 
as  the  measurement  of  the  nuclear  quadrupole  moment  Q)  is  complicated  by  the 
contributions  to  the  hyperfine  interactions  arising  from  the  distortion  of 
the  otherwise  spherical  closed  electronic  shells  of  the  system.  One  of 
Sternheimer ' s  important  contributions  was  the  striking  demonstration  that, 
for  an  ion  having  a  nuclear  quadrupole  moment,  the  quadrupole  interaction 
arising  from  the  field  induced  (1)  by  external  charges  (as  in  a  salt),  or 
(2)  from  the  ion’s  own  aspherical  charge  distribution  (if  the  ion  is  not 
spherically  symmetrical),  was  changed  appreciably  by  the  distortion  of  the 
ion’s  closed  shells.  Following  Sternheimer ' s  pioneering  investigation1, 
much  work  has  been  done  using  the  perturbation-variation  techniques  which 
he  developed;  more  recently  Das  and  Bersohn  have  used2  a  (less  precise) 
analytic  function,  instead  of  the  more  laborious  numerical  integration 
procedure  in  solving  the  perturbation  theory  equations. 

In  this  paper3,  we  develop  a  method  for  calculating,  within  the  Hartree- 
Fock  formalism,  the  distortions  induced  in  electronic  distributions  by 
valence  electrons  and/or  crystalline  fields  and  their  effects  on  magnetic 
and  electric  hyperfine  interactions.  For  illustrative  purposes  emphasis 
is  placed  on  the  calculation  of  Sternheimer  quadrupole  polarizabilities  and 
anti-shielding  factors,  and  some  results  are  reported  for  the  anti-shielding 
factors  (y^)  of  Cu+  and  Cl*.  Since  a  self-consistent  field  procedure  is 
followed,  the  distortions  induced  in  the  inner  closed  shells  by  the  distorted 
outer  closed  shells  are  included  in  a  natural  way  and  by  comparison  with 
the  results  of  the  perturbation  method  (which  does  not  take  these  into 
consideration)  these  additional  effects  are  shown  to  be  significant.  The 
problems  of  orthogonality,  exchange,  and  self-consistency  which  have 
complicated  applications  of  the  perturbation  method  are  easily  resolved  by 

this  approach.  Although  the  present  results  are  for  spherical  ions  in  an 

external  field  our  method  is  also  applicable  to  the  problem  of  induced 

electric  quadrupole  (and  magnetic  dipole  and  higher  multipole)  distortions 

of  an  ion  by  its  own  aspherical  electronic  charge  distributions  (e.g. ,  3d 
or  4f  electrons),  a  problem  recently  highlighted  by  MSssbauer  measurements. 

In  contrast  to  the  present  approach,  perturbation  theory  calculations  for 
these  are  difficult  to  carry  out. 


II.  PERTURBATION-VARIATION  METHOD 

The  perturbation  theory  approach  to  the  problem  of  the  polarization  of 
electron  shells  of  atoms  and  ions  has  been  discussed  previously  at  length1’*'3. 
Some  brief  details  will  be  given  here  in  order  to  provide  a  framework  for 
discussing  our  method  and  for  making  comparisons 
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Writing  the  Hamiltonian  in  the  uaual  form,  H  *  H0  +  H|,  with  He  the 
unperturbed  Hamiltonian  and  Hx  a  perturbing  potential,  the  first  order 
perturbation  i (r )  to  the  unperturbed  solution,  \ p0  (r),  is  determined  from 
the  relation 

(H0  -  E0)^(r)  -  -  (Hi  -  Ejl^jfr)  (1) 

which  was  obtained  by  applying  the  variational  principle  to  the  perturbation 
equations.  Here 

E0  “  <'/'„IhJ’/'0>.  Ex  =  <00 1  Hx  l'/'0> .  and  >  -  0. 

Equation  1  has  been  solved  in  two  distinct  ways:  (a)  by  exact  numerical 
solution,  as  was  done  by  Sternheimer  and  collaborators1,  and  (b)  by, the 
analytic  approach  of  Das  and  Bersohn2  in  which  the  radial  part  of  ip\(r)  is 
assumed  to  be  related  to  the  radial  part  of  0°(r)  by 

<(r)  ■  u'0(r)  §  aB  r"  (2) 

and  the  parameters  aa  are  determined  by  minimizing  the  energy  with  respect 
to  variation  of  these  parameters.  This1  technique  has  the  advantage  of  being 
easier  to  carry  out  than  Sternheimer ’ s  but  is  inferior  because  full  varia¬ 
tional  freedom  is  not  accorded  to  u'x  (e.g.,  if  uB  is  a  noded  function  then 
u^  is  constrained  to  have  the  same  nodes).  In  more  recent  applications  the 
problems  of  properly  maintaining  orthogonality6,  self-consistency7,  and  in¬ 
cluding  exchange8,  have  been  emphasized. 

Since  for  quadrupole  interactions  the  perturbing  potential  Hx  has 
(6 ,4>)  symmetry,  perturbed  orbital  character  will  be  mixed  into  the  unper¬ 
turbed  orbitals,  i/*8,  in  the  following  ways: 

0J-N.  .O/'J  ♦  ’/'j ) 

0p-NpC0j  +  0p  +  0f  >  4  (3) 

.  2  ♦  0Je  ♦  2);  ^ 1 

where  the  Nx  are  normalization  constants.  The  mixing  of  a  41'  of  -t  in  common 
with  i is  called  a  "radial"  excitation  whereas  the  (/>'  components  having 
't  i  2  are  called  "angular"  excitations. 

Sternheimer  first  discussed  the  two  cases  already  referred  to:  (1)  Hx 
arising  external  to  the  ion  and  (2)  H,  due  to  an  aspherical  unclosed 
electron  shell  within  the  ion.  Case  (2)  and  in  particular  the  associated 
magnetic  hyperfine  effects  will  be  discussed  in  subsequent  papers;  in  what 
follows  we  shall  outline  our  method  an<j  present  results  only  for  case  (1). 
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For  a  quarupole  potential  due  to  charges  external  to  the  central  ion, 
the  interaction  with  the  nuclear  quadrupole  moment,  Q,  ia  given  aa  e2  q  Q, 
where  q  ia  the  external  charge  electric  field  gradient  at  the  nucleus, 
Sternheimer  found  that  the  external  potential  distorted  the  electron  shells 
of  the  ion  and  that  this  distortion  also  interacted  with  Q  giving  a  total 
interaction  of  the  form 


HQ  -  e2  q  Q  (1  •  yj  (4) 

Here  ym  is  the  Sternheimer  antishielding  factor.  For  ions  consisting  of 
more  than  closed  s  shells,  y9  is  negative  and  large  (typicadlly  *10  to  >100), 
The  "radial*  contributions  to  ym  are  generally  much  larger  than  (and  of 
opposite  sign  to)  the  "angular"  contributions. 


III.  HARTREE-FOCK  SELF-CONSISTENT  FIELD  METHOD 

AND 

QUAORUPOLE  POLARIZABILITIES 

Spin  (or  exchange)  polarized  Hartree-Fock  calculations  have  recently 
played  a  prominent  role  in  understanding  magnetic  hyperfine  interactions9. 

As  has  been  previously  discussed,  such  calculations  involve  relaxing  th>e 
restriction  imposed  on  conventional  atomic  Hartree-Fock  solutions  that 
electrons  in  the  same  shell  but  differing  in  spin  (i.e.,  m,  quantum  number) 
have  the  same  radial  wave  function. 

The  mB  restriction  is  not  the  only  one  imposed  on  H-F  solutions9.  Two 
other  constraints,  which  concern  us  here,  are:  (A)  that  the  spatial  part 
of  the  one-electron  functions  be  separable  into  a  radial  function  times  and 
angular  function  (for  free  ions  this  is  a  single  spherical  harmonic)  and  (B) 
that  electrons  of  the  same  shell  but  of  differing  magnetic  quantum  number 
(m^)  have  the  same  radial  wave  function.  These  restrictions  lead  to  the 
shell  structure  description  of  atoms,  molecules,  and  solids  but  are  only 
rigorously  valid  for  closed  shell  atoms;  their  imposition  gives  not  only  a 
physically  simple  picture  of  electronic  systems  but  also  allow  the  computa¬ 
tions  to  be  of  reasonable  magnitude.  By  analogy  with  the  well-known  spin- 
polarized  H-F  calculations  one  may  conveniently  define  calculations  in 
which  restriction  (B)  is  relaxed  as  "orbitally"  polarized  H-F  calculations. 

We  have  investigated  the  physical  consequences  of  relaxing  conditions 
(A)  and  (B)  and  have  found  that  quadrupole  polarizabilities  and  Sternheimer 
anti -shielding  factors  can  be  obtained  by  these  means.  Relaxing  (A)  results 
in  "angular"  distortions  of  the  electron  shells  whereas  relaxing  (B)  yields 
"radial"  distortions.  (Gq.  3  was  written  so  as  to  make  this  appear  more 
self-eviddnt» /.  The  method  is  applicable  to  the  distortions  due  to  either 
the  case  of  an  external  perturbing  field19 or  that  arising  from  an  ion’s 
aspher  ical  charge  distribution.  For  a  purely  ionic  crystal  "potential", 
such  as  the  one  we  discuss  in  this  paper,  only  electric  hyperfine  effects 
are  induced11;  the  case  of  induced  magnetic  and  electric  hyperfine  ineractions 
due  to  an  unclosed  shell  with  the  ion  will  be  discussed  in  a  future  paper. 
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A.  Method  of  Computation 


For  simplicity,  consider  a  closed  shell  ion  (e.g.,  Cl*  and  Cu+)  in  the 
field  of  a  single  external  point  charge  a  distance  R  away.  In  atomic  units, 
the  quadrupole  term  of  the  potential  expanded  in  spherical  harmonics  and 
powers  of  r  (the  distance  from  the  nucleus  having  quadrupole  moment,  Q,  to 
some,  point  P),  is  simply 


VQ  -  -  Zv/I75  ®  (2,0)  ll  (5) 

R3 

where  Z  is  the  magnitude  of  the  point  charge  and  0  (2,0)  is  the  normalized 
spherical  harmonic  defined  by  Condon  and  Shortley.  In  otder  to  carry  out 
the  computations,  we  include  Vq  into  our  free  ion  Hamiltonian  and  rederive 
the  Hartree-Fock  equations.  These  equations  are  then  solved  (in  principle) 
in  a  straightforward  way;  in  practice,  only  the  radial  anti-shielding  terms 
have  been  obtained  as  this  can  be  done  with  existing  H-F  computational 
machinery.  The  less  important  angular  distortions  require  a  H-F  treatment 
involving  functions  of  mixed  angular  character  and  can  be  obtained  if  one 
utilizes  and  extends  existing  analytic  H-F  techniques12.  Since  any  errors 
associated  with  crude  estimates  of  angular  anti -shie lding  tend  to  be  dwarfed 
by  inadequacies  associated  with  the  crystal  field  model  used  in  ym  investiga¬ 
tions  (to  be  discussed  later),  the  effort  necessary  to  extend  the  H-F 
computational  machinery  to  include  angular  anti -shielding  did  not  seem 
warranted  here. 

The  computations  were  done  with  analytic  H-F  methods  used  and  described 
previously13'14.  A  quadrupole  field  due  to  an  external  charge  distribution 
was  included  in  the  Hamiltonian  from  which  the  H-F  equations  were  derived, 
subject  to  restriction  (B)  being  relaxed.  In  the  solution  of  these  equa¬ 
tions,  the  ion's  shells  were  self-consistently  distorted  and  the  quadrupole 
interaction  of  these  distortions  with  the  ion's  nucleus  then  yielded15^. 

Using  this  scheme,  calculations  have  been  done  for  Cl*  and  Cu+.  The 
anti -shie lding  for  these  ions  is  of  considerable  interest;  numerical  and 
analytic  perturbation  results  already  exist  so  that  comparisons  may  be 
readily  made.  In  the  course  of  the  investigation,  a  conventional  H-F  cal¬ 
culation  was  also  done  for  Cu+  which  is  of  greater  accuracy  than  those 
appearing  in  the  literature.  The  resulting  one-electron  functions  are 
tabulated  in  Appendix  I.  The  basis  set  appearing  in  this  calculation  was 
used  in  the  ym  calculation.  A  previously  published  conventional  H-F14 basis 
set  was  used  for  Cl*. 

B.  Results 

The  radial  contributions  to  ym  as  given  by  SCF  and  perturbation  theory 
calculations  for  Cu+  and  Cl*  are  listed  in  Table  I.  The  three  sets  of 
results  are  seen  to  be  in  very  good  agreement  for  Cu+  but  differ  markedly 
for  Cl*.  Perhaps  the: most  significant  feature  of  the  results  is  the  increase 
in  y 2 p_p  by  a  factor  of  two  for  Cu+  and  eight  for  Cl".  These  increases 
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arise  from  the  distortions  produced  in  the  inner  2p  shells  by  the  distorted 
outer  3p  (and  3d)  shells  -  distortions  which  are  included  in  our  SCF  treat¬ 
ment  but  not  in  the  perturbation  theory  approach.  Upon  reflection  it  may 
appear  somewhat  surprising  that  the  inner  shell  enhancement  is  not  even 
grester  for  this  shell  in  view  of  the  large  distortions  of  the  outer  shells 
produced  by  the  applied  field.  However,  the  are  not  enhanced  by 

factors  proportional  to  the  y3  ^p‘s  because  unlike  the  nucleus,  the  2p 
shell  is  overlapped  by  the  3p  (and  3d)  shell.  As  soon  as  the  perturbing 
charge  overlaps  the  perturbed  shell  one  must  replace  the  point  charge 
potential  given  by  Eq.  5  by 

V  (r)~/  p(x)  Jil  dx  ♦  fp(x)  JlL  dx  (6) 

’  •  r3  r  x3 

where  p(x )  is  the  radial  distribution  of  the  perturbing  charge  density. 

Such  a  distributed  charge  yields  a  smaller  y2p_  than  would  either  a  Vq 
given  by  Eq.  5  or  one  arising  from  a  distributed  charge  completely  external 
to  the  2p  shell.  In  fact,  a  perturbation  theory  estimate  of  the  Cl*y2  , 
utilizing  Eq.  6  and  a  p(x)  consisting  of  the  external  charge  plus  the  dis¬ 
torted  3p  shell,  crudely  reproduces  (to  better  than  a  factor  of  two)  the 
value  appearing  in  Table  I.  In  other  words,  where  anti -shielding  is  con¬ 
cerned,  a  shell  such  as  the  2p  is  far  from  being  completely  *inside*  a  3p 
(or  3d)  shell.  (If  it  were  otherwise  we  would  have  a  Cl”  72 p-p  ~‘100.) 


TABLE  I 

COMPARISON  OF  "RADIAL"  CONTRIBUTIONS  TO  ym 


Second  Order  Perturbation  Theory 

Cl* 

Present  Results 

Numerical* 

Analytical- 

variational** 

2P 

-12.1 

-  1.5 

~  -1 

3p 

-78.3 

-56.5 

-50  to  -60 

Cu+ 

2p 

-  1.2 

-  .62 

-  .6 

3P 

-  7.8 

-7.9 

-5.8 

3d 

-  8.6 

-8.5 

-8.6 

If  one  assumes  that  the  enhanced  SCF  7^,%  values  are  entirely 
associated  with  self-consistency  (i.e.,  an  aspherical  potential  which  in¬ 
cludes  the  effects  of  the  closed  shell  distortions),  one  can  use  the 
results  of  Table  I  to  estimate  the  same  effect  on  the  angular  7  terms. 
Results  are  given16  in  Table  II  where  estimates  of  the  total  ym's  are  also 
listed.  One  could  instead  directly  estimate  the  angular  Sternheimer  anti¬ 
shielding  utilizing  Eq.  6;  for  Cl”  this  method  gave  the  ya  (angular)  value 
of  3.4  appearing  in  Table  II. 


*R.M.  Sternheiner  and  B.M.  Poley,  Phys.  Rev.  IQ£,  731  (1956). 
,,t.O.  Vikner  and  f.P.  Das,  Phys.  Rev.  109,  380  (1958). 


-7- 


TABLE  II 

RADIAL  AND  ANGULAR  CONTRIBUTIONS  TO  FOR  Cl'  AND  Cu+ 


cr 

Cu+ 

Present 

Results 

Second  Or dor 
Perturbation 
Theory 

Present 

Results 

Second  Order 
Perturbation 
Theory 

ym  (radial) 

•90.4 

•58.0 

-17.8 

-17.0 

(angular) 

3.4  to  5.5 

1.5* 

~  0.6 

~  0.1 

TOTAL  y^ 

-87.0  to  -84.9 

•56.5 

-17.0 

-16.9 

It  is,  of  course,  naive  to  expect  that  the  differences  between  the 
perturbation  theory  and  SCF  results  of  Table  I  are  due  entirely  to  self- 
consistency.  Khubchandani ,  Sharma ,  and  Das17,  using  a  modification  of 
Dalgarno's  method5  for  including  exchange,  showed  that  the  omission  of 
exchange  in  the  standard  perturbation  calculations  can  be  numerically 
significant  and  Ingall’s  has  shown4  that  the  ma intainance  of  orthogonality 
can  have  severe  effects  og  a  computed  ym .  Let  us  now  examine  the  most 
noticeable  feature  of  the  results  listed  in  Table  I,  i.e.,  the  substantially 
larger  value  for  the  Cl"  3p  terai  ^3p-p)  *«  have  obtained.  In  order 

to  determine  the  cause  of  this  difference  we  have  carried  out  a  series  of 
analytic  perturbation  calculations  for  an{*  ^3p-p  malntaln^n8 

orthogonality  in  a  variety  of  ways.  From  the  results  we  conclude  that  the 
increased  SCF  value  for  73p^p  is  largely  due  to  the  proper  maintainance  of 
orthogonality  with  the  2p  shell.  (These  computations  also  agree  with 
Ingall's  observation  that  the  results  are  extremely  sensitive  to  the  way 
that  orthogonalization  is  carried  out;  details  will  be  given  in  a  future 
publication.  ) 

A  detailed  breakdown  of  the  respective  roles  of  self -consistency, 
exchange  and  orthogonality  has  not  been  attempted  but  it  appears  that  the 
former  is  important  for  inner  shells  and  that  the  latter  contributes 
heavily  to  the  final  character  of  outer  shell  terms. 


IV.  DISCUSSION 

Orthogonality,  self-consistency,  and  exchange  are  naturally  and 
properly  handled  by  the  method  we  are  proposing;  in  addition,  the  method 
is  computationally  practicable.  For  these  reasons,  we  believe  its  use 
offers  important  advantages  over  the  traditional  methods  for  obtaining  ym 1 s . 

A  major  defect  of  any  of  the  methods  for  obtaining  y m  lies  in  the 
deficiencies  introduced  by  the  approximate  form  of  the  assumed  perturbing 
potential.  For  a  given  external  potential  for  which  a  ym  has  been  computed, 
there  still  remains  the  question  of  agreement  with  a  y  derived  somehow  from 
qQ(l  -  y)  gotten  from  experiment.  The  absence  of  exact  values  of  q  and  Q 


•Jf.JV.  Sternheiner  (unpublished) . 
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make  it  impossible  to  make  definite  statements  concerning  this  question 
but  it  has  been  suggested16  that  y'  s  of  the  order  of  -10  to  -15  are  appro¬ 
priate  for  Cl*  and  Cu+  implying  that  the  calculated  value  of  ym  for  Cl*  is 
a  severe  overestimate  whereas  the  Cu  +  value  is  only  slightly  so.  This 
comparison  raises  doubts  concerning  the  accuracy  of  a  computed  q  and  the 
related  V  ,  which  go  beyond  the  question  of  defining  a  potential  in  that 
it  involves  the  concept  of  an  "ion”  in  a  solid  and  of  the  "potential" 
associated  with  it.  It  may  be  argued  that  such  a  model  is  inadequate  for 
the  problems  of  interest  to  us  here;  let  us  consider  this  matter. 

The  lattice  sum  over  ion  point  charges,  which  is  normally  done  when 
estimating  an  electric  field  gradient,  yields  a  potential  of  the  form  of 
Eq.  5.  Although  the  results  may  be  adequate 19  for  obtaining  the  direct 
quadrupole  interaction  with  the  nucleus,  such  a  sum  leaves  much  to  be 
desired  when  used  as  the  perturbing  potential  seen  by  the  ion's  electrons 
because  the  ion's  charge  density  overlaps  that  of  the  neighboring  ion 
electrons  and  nuclei.  Hence,  a  potential  of  the  form  r2Y2(0,$)  is  clearly 
not  appropriate  over  that  region.  In  addition  to  more  accurately  defining 
an  electrostatic  field,  problems  associated  with  wave  function  orthogonality 
(between  neighboring  ions)  and  covalency  (if  some  of  the  ions  have  unclosed 
shells)  occur.  Also  the  crystal  potential  will  have  components  of  other 
than  quadrupole  symmetry  which  may  distort  the  ion  and  in  turn  affect  the 
electronic  quadrupole  distortion.  For  example,  it  has  been  argued  that  the 
tendency  of  negative  ions  to  contract  when  inserted  into  ionic  crystals  may 
possibly  result  in  severe  repercussions18  for  y_*. 

The  importance  of  some  of  these  problems  can  be  illustrated  if  we 
consider  y  computed  as  a  function  of  the  distance  (R)  of  the  source  of  q 
from  the  nucleus  in  question;  utilizing  Eq.  6,  with  p(x)  •  8(x  -  R),  one 
obtains  Sternheimer ’ s  7(R)  factor1.  The  y3p_p  result  for  Cl  is  plotted 
as  the  ratio  y  (R) /yn  in  Figure  1.  Since  internuclear  spacing  are  typically 
of  the  order  of  5  a.u.,  the  arguments  of  the  previous  paragraph  are  made 
more  apparent  by  the  variation  of  y(R)  with  R  which  is  seen  in  the  figure. 

At  best,  this  means  that  nonquadrupole  crystal  field  effects  and  ortho- 
gonalization  repercussions  must  be  dealt  with  before  a  realistic  value  is 
obtained  for  y.  At  worst,  it  may  mean  that  the  tight  binding  approach  of 
perturbing  an  essentially  free  ion  by  a  crystalline  environment  must  be 
abandoned  before  a  substantial  improvement  over  current  yw  estimates  is 
made 
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APPENDIX  I 


A  conventional  H-F  Cu+  wave  function,  superior  to  those  in  the 
literature,  is  tabulated  here  whose  accuracy  is  similar  to  a  previously 
reported2flMn2+  wave  function.  Orthonormal  analytic  Hartree-Fock  orbitals, 
Ui  (r),  of  the  form 

Ui  <r>  ’  Sj  Cij  Rj  <r>  A(l) 

were  obtained,  normalized  such  that 


J  |Ui  <r)  | 2  dr  -  1. 

O 


A  ( 2  > 


The  basis  functions,  Bj ,  are  of  the  form 


Rj  (r) 


N ; 


(l  *  A;  +1> 


-2jr 


A  ( 3 ) 


where  Nj  is  a  normalization  constant  and  is  expressible  in  terms  of  the 
other  parameters,  i.e., 


a  +  3 

[(#j  )  +  2Aj  /ill  *  2Aj  ♦  2)!] 


1/2 


A  ( 4 ) 


Ujlr)  of  common  l  value  are  constructed  from  a  common  set  of  R.(r)’s. 

Z  j  *  s ,  AVs  and  combining  coefficients  (C^.’s)  are  listed  in  Table  III;  the 
one-electron  energies  are  given  in  Table  iV.  The  C^'s  define  orthonormal 
functions  to  the  number  of  digits  reported  but  have  by  no  means  been 


uniquely  determined  to  that  number  of  digits.  The  total  energy  for  the  ion 
is  -1638.724.  a.u.  as  compared  with  -1638.705  a.u.  obtained  in  a  previous 
calculation  .  The  improvement  is  primarily  associated  with  the  greater 
variational  freedom  allowed  to  the  outer  parts  of  the  3s  and  3p  shells; 
repercussions  on  the  3d  shell  were  small. 
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TABLE  III 


PARAMETERS  (AB  AND  Za)  AND  COMBINING  COEFFICIENTS 
(Cjj)  DEFINING  THE  HARTREE-FOCK  ORBITALS  FOR  Cu+ 


j 

aj 

zi 

C1..J 

C2i.j 

C3«.j 

1 

0 

84.8828 

.00037044 

-.00036468 

-.00038449 

2 

0 

30.4491 

.91205876 

-.27950925 

.10652016 

3 

1 

26.9961 

.10841891 

-.16245115 

.04976048 

4 

1 

14.1023 

-.00647884 

.64460027 

-.21227812 

5 

2 

13.3297 

.00561417 

.46415723 

-.38260515 

6 

2 

8.5681 

-.00204424 

.06610429 

.17263958 

7 

2 

5.0982 

.00065755 

-.004  54834 

.94245212 

8 

2 

2.8171 

-.00016921 

.00134453 

.07S47339 

C2p.  j 

C3p.j 

0 

34.S244 

.00103326 

.00096650 

0 

18.9852, 

.13481246 

-.05841121 

11 

0 

11.6372 

.84951100 

-.29915630 

12 

1 

10.8502 

.02041835 

-.11625686 

13 

1 

6.8083 

.02746109 

.54128259 

14 

1 

3.9668 

-.00525461 

.63298421 

15 

1 

2.2563 

.00172020 

.02228438 

C3d.  j 

16 

a 

1.6956 

.19657400 

| 

n 

2.4954 

.25771318 

u 

4.2881 

.46579549 

7.5015 

.24455545 

53 

El 

13.0768 

.02287090 

TABLE  IV 

ONE -ELECTRON  ENERGIES  FOR  Cu+  IN  ATOMIC  UNITS 
(1  a.u.  ■  2  ry) 
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